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摘 要 : 陕西 省 消 北 矿区 地 处 黄土 高 原 ,受到 自然 环境 和 煤炭 开采 等 人 类 活动 的 双重 影响 ,该 地 
土壤 侵蚀 更 为 严重 ,生态 环境 更 加 脆弱 。 以 渭 北 矿 
壤 侵 人 蚀 进 行 了 评估 ,并 综合 植被 覆盖 度 、 多 年 平均 降水 量 、 坡 度 、 土 地 利用 类 型 和 煤炭 年 产量 等 影响 
妹子 ,应 用 地 理 探 测 器 方法 对 渭 北 矿区 土壤 侵 保 进 行 定量 归 因 。 研 究 结果 对 矿区 水 土 流 失 防 治 其 
有 一 定 参 考 价 值 。 结 果 表 明 :(1) 涓 北 矿区 土壤 侵蚀 以 微 度 和 轻 度 侵蚀 为 主 , 土 壤 侵蚀 严重 的 区 域 


FFA, E 3X 
徐州 ”221116) 


I4 


区 为 研究 对 象 ,基于 RUSLE 模型 对 渭 北 矿区 土 


> 


主要 位 于 研究 区 西南 部 、 中 部 和 东南 部 。(2) 植 被 覆盖 度 和 多 年 平均 降水 量 是 造成 研究 区 土壤 侵 
蚀 的 主导 因子 ,坡度 介 于 20° ~25° 的 地 区 、 植 被 覆盖 度 小 于 0.3 的 区 域 和 裸 地 是 发 生 土壤 侵蚀 的 高 
风险 区 。(3) 渭 北 矿区 各 因子 协同 作用 对 土壤 侵蚀 的 解释 力 均 大 于 单 因子 解释 力 , 因 此 多 个 因子 


共同 作用 会 对 土壤 侵蚀 造 成 显著 影响 。 
关键 词 : 


位 于 陕 北 黄土 高 原 区 南 缘 的 渭 北 矿区 是 我 国 重 


土壤 侵蚀 ; RUSLE 模型 ; 黄土 高 原 ; 地 理 探 测 器 ; 煤矿 区 ; 驱动 因素 


区 域 进行 了 预测 ; 秦 伟 等 ”利用 RUSLE 模型 评估 了 


要 的 煤 痰 生产 基地 。 有 上 脆弱 的 生态 环境 和 高 强度 的 人 
类 活动 使 得 该 地 区 的 土壤 侵蚀 现象 日 益 加 剧 ,对 生 
态 环境 造成 了 严重 破坏 ,成 为 社会 经 济 发 展 的 重 
要 障碍 。 因 此 ,定量 评价 该 区 域 的 土壤 侵蚀 现状 , 科 
学 认识 和 定量 评估 地 理 环 境 因 子 `、 人 类 活动 及 其 耦 


陕西 省 四 面 窑 沟 流域 的 土壤 侵蚀 强度 ,所 得 结果 与 
实际 调查 数据 较为 吻合 ,说 明基 于 RUSLE 模型 评估 
黄土 区 小 流域 土壤 侵蚀 的 方法 有 效 、 可 行 。 在 对 区 
域 土壤 侵蚀 进行 定量 分 析 的 基础 上 进一步 探究 其 驱 
动因 素 ,能 够 帮助 明确 水 土 流失 治理 方向 。 在 土壤 


合作 用 对 土壤 侵蚀 的 影响 ,对 矿区 水 土 保持 和 环境 
保护 具有 重要 意义 。 

与 传统 的 土壤 侵蚀 调查 方法 相 比 ,利用 RS 和 
GIS 技术 进行 相关 地 理 信 息 的 搜集 及 分 析 , 不 仅 提 
高 了 工作 效率 ,更 使 得 科学 性 和 准确 性 有 了 保 
障 '"。 在 大 量 土壤 侵蚀 试验 和 定量 观测 数据 的 基 
础 之 上 ,通过 不 断 完善 改进 形成 的 修正 的 通用 土壤 
流失 方程 (Revised Universal Soil Loss Equation , RU- 
SLE) 是 目前 在 土壤 侵蚀 定量 估算 方面 应 用 最 广泛 
的 经 验 坡 面 土壤 侵蚀 预报 模型 ”。 近 年 来 ,该 模型 
已 经 被 成 熟地 运用 到 我 国 黄土 高 原 地 区 的 土壤 侵蚀 
研究 当中 。 李 柏 延 等 中 基于 RUSLE 模型 对 榆林 市 
2001—2010 年 连续 10 a 土壤 侵蚀 动态 变化 及 其 分 
布 进行 了 分 析 ,并 对 未 来 几 年 土壤 侵蚀 发 生 的 主要 
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侵蚀 驱动 因子 识别 方面 , 周 来 等 "利用 闪 加 统计 分 
析 和 相关 性 分 析 的 方法 ,探究 了 福建 省 将 乐 县 坡度 
和 土地 覆盖 类 型 与 土壤 侵蚀 强度 之 间 的 相关 关系 ， 
指出 15° ~25° 坡 度 带 的 土壤 侵蚀 量 占 比 最 大 并 且 
森林 和 稀 玻 植被 对 土壤 保持 有 着 较 大 的 贡献 ; 王 
涛 ”以 自然 状态 下 陕 北 洛 河流 域 土壤 侵蚀 模 数 . 退 
耕 还 林 还 章 工 程 实施 区 域 自 然 状 态 土壤 侵蚀 模 数 等 
数据 为 基础 ,采用 礁 加 统计 分 析 的 方法 ,定量 评估 了 
退耕 还 林 还 草 工程 的 实施 和 降雨 对 土壤 侵蚀 的 影响 。 

通过 对 文献 的 梳理 可 以 发 现 ,目前 对 于 黄土 高 
原 区 的 矿区 土壤 侵蚀 定量 评价 研究 还 较 少 。 作 为 陕 
西 省 开采 历史 悠久 的 重要 煤炭 基地 ,反映 渭 北 矿区 
土壤 侵蚀 现状 的 研究 还 鲜 有 报道 。 此 外 , 现 有 的 工 
作对 于 土壤 侵蚀 驱动 因素 的 分 析 虽 然 均 聚 焦 于 自然 
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环境 和 人 类 活动 的 影响 ,但 多 采用 相关 性 分 析 或 回 
归 统 计 等 常规 方法 ,而 对 多 因子 间 交 互 作用 进行 定 
量 分 析 , 衡 量 其 对 土壤 侵蚀 影响 力 的 研究 还 较 少 。 
基于 地 理 探 测 器 的 自身 优势 ,利用 该 方法 可 以 通过 
探测 地 理 现象 的 空间 分 层 异 质 性 来 揭示 其 背后 的 驱 
动力 ,能够 有 效 地 识别 地 理 现 象 的 高 风险 区 以 及 各 
风险 因子 间 的 交互 作用 。 因 此 本 文 在 RS 和 GIS dx 
术 的 支持 下 ,基于 RUSLE 模型 对 渭 北 矿区 的 土壤 侵 
蚀 现 状 进行 定量 评价 并 利用 地 理 探测 器 的 方法 开展 
渭 北 矿区 土壤 侵蚀 多 因子 定量 归 因 ,以 期 甄别 该 区 
域 土壤 侵蚀 的 主导 影响 因子 ,分 析 各 因子 间 的 协同 
作用 ,为 矿区 水 土 保 持 治理 和 环境 保护 政策 提供 科 


1 研究 区 概况 


渭 北 矿区 位 于 陕西 省 中 部 渭河 以 北 ,由 西向 东 
主要 分 为 铜川 蒲 白 、 溢 合 、 韩 城 4 个 矿区 ,是 陕西 省 
主要 产 煤 区 之 一 ,地 理 位 置 为 34°38’ ~35°47'N, 
108*40' ~110°37'E, 面 积 约 8 000 km , 占 全 省 面积 
的 4.86% (图 1)。 矿区 地 貌 类 型 兼 有 平原 和 丘陵 
的 特点 ,地貌 较为 复杂 ,小 起 伏 中 低 山 \ 中 低 海拔 台 
地 、 低 高 海拔 平原 .丘陵 兼 而 有 之 。 矿 区 地 形 呈 现 南 
低 北 高 的 特点 ,坡度 起 伏 较 大 ,最 小 坡度 为 0° ,最 大 
坡度 为 70"。 矿 区 内 土壤 大 部 分 颗粒 组 成 是 粉 沙 粒 


径 ,质地 均 细 ,组 织 疏 松 , 土 粒 间 主要 靠 碳酸 盐 胶 结 ， 
极 易 在 水 中 前 解 分 散 , 抗 人 蚀 力 薄弱 。 黄 河 最 大 的 
支流 一 一 渭河 位 于 矿区 南部 ,矿区 内 部 河流 支流 众 
多 ,对 土壤 的 冲刷 作用 造成 区 域内 沟 蜜 纵横, 地形 破 
碎 , 沟 秘密 度 普遍 达到 1.0 km * km 左右。 矿区 内 
矿藏 储量 丰富 ,煤炭 资源 量 为 371. 12 x 10* 1t, 约 占 
陕西 省 煤炭 资源 总 量 的 10% , 煤 种 齐全 , 煤 质 优良 ， 
是 陕西 省 最 主要 的 产 煤 区 。 特 殊 的 地 势 .土壤 类 型 
等 自然 条 件 和 强烈 的 人 类 活动 使 得 渭 北 矿区 土壤 侵 
蚀 状况 不 容 乐观 。 根 据 我 国 水 利 部 公布 的 土壤 侵蚀 
分 类 分 级 标准 ,该 地 区 土壤 侵蚀 类 型 以 水 力 侵蚀 为 
主 ,侵蚀 强度 主要 为 微 度 和 轻 度 侵蚀 。 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 及 处 理 

基于 RUSLE 模型 ,对 渭 北 矿区 土壤 侵蚀 现状 进 
行 评估 ,所 需 人 研究 数据 如 下 :土壤 数据 主要 来 源 于 两 
个 方面 ,一 是 由 中 国 科学 院 资 源 环 境 科学 数据 中 心 
(http: /www. resdc. cn) 提供 的 全 国 范围 土壤 砂粒 、 
粉 粒 和 黏 粒 的 百分比 含量 ,二 是 由 国土 部 门 提 供 的 
土壤 采样 点 数据 和 土壤 类 型 数据 ,并 通过 插值 生成 
研究 区 土壤 有 机 碳 分 布 图 ;地 形 数据 以 来 源 于 中 国 
科学 院 计算 机 网 络 信息 中 心 国际 科学 数据 镜像 网 站 
( http : /www. gscloud. cn ) 的 GDEMDEM 30 m 分 辨 
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1 渭 北 矿区 地 理 位 置 示意 图 


Fig.1 Geographical location of the Weibei mining area 
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率 数字 高 程 数据 为 基础 , 拼接 裁剪 生成 研究 区 
DEM ,并 利用 ArcGISIO. 2 软件 中 表面 分 析 模 块 提取 
坡度 等 数据 ;多 年 平均 降水 数据 和 月 值 降 水 数据 来 
源 于 中 国 气象 数据 网 (http: // data. cma. en) ,选取 渭 
北 矿区 及 其 周边 区 域 31 个 气象 站 点 的 降水 量 数据 
插值 生成 覆盖 研究 区 的 多 年 平均 降水 量 栅 格 数据 ; 
植被 覆盖 数据 来 源 于 美国 国家 航空 航天 局 (NASA ) 
官网 (http: // ladsweb. nascom. nasa. gov ) ,采用 2018 
4E 6 -9 月 MODIS 数据 250 m 分 辨 率 16 d 合成 产品 
MODI3Q1 影像 ,利用 ENVI 5. 1 软件 计算 矿区 归 一 
化 植被 指数 (NDVI) 值 ,并 利用 像 元 二 分 模型 计算 植 
被 覆盖 度 ,其 值 范 围 在 0 ~ 1; 土 地 利用 数据 由 2018 
年 3 月 和 4 月 的 四 景 Landsat-8 OLI 影像 监督 分 类 
得 到 ,利用 ENVIS. 1 软件 进行 人 机 交互 解 译 得 到 
30 m 分 辨 紊 土地 利用 类 型 数据 并 进行 精度 验证 ,总 
体 分 类 精度 为 96.58% ,数据 较为 可 靠 ;煤矿 数据 由 
2017 年 和 2018 年 实地 调研 资料 整理 得 到 ,包括 滑 
北 矿 区 矿 点 名 称 、 生 产 状况 和 煤炭 年 产量 等 ,用 于 反 
映 矿区 的 开采 活动 。 

将 上 述 各 类 栅 格 数据 统一 为 30 m x30 m 像 元 
大 小 ,并 统一 采用 WGS_1984_Lambert_Conformal 
Conic 投影 。 

2.2 研究 方法 
2.2.1 RUSLE 模型 RUSLE 模型 是 在 美国 通用 
土壤 流失 模型 的 基础 上 发 展 得 来 的 ,该 模型 是 目前 
国内 外 应 用 广泛 的 土壤 侵蚀 预测 模型 之 一 。RUSLE 
模型 不 仅 弥 补 了 实地 观测 在 大 尺度 应 用 上 的 局 限 
性 ,而 且 在 不 同 尺度 的 模拟 中 取得 了 较 好 的 效 
RU. py, 4t RUSLE 模型 作为 研究 渭 北 矿区 
土壤 侵蚀 量 的 定量 模型 ,其 表达 式 为 : 
A-RxKxLxSxCxP 


(1) 


式 中 :4 为 年 土壤 侵蚀 量 (t' hm 7 < a7) ;R 为 降雨 
侵蚀 力 因 子 (MJ - mm. hm? h^ +a); KHE 
壤 可 蚀 性 因子 (t + hm? + h + MJ + mm™' + hn ); 
了 为 坡 长 因子 ;$ 为 坡度 因子 (°*);C 为 地 表 植 被 覆 
HAF; P 为 水 土 保 持 措施 因子 。 

降雨 侵蚀 力 是 指 由 于 降雨 引起 的 土壤 潜在 的 侵 
蚀 能 力 ,是 造成 区 域 土壤 侵蚀 的 主要 因素 之 一 。 根 
据 数据 获取 情况 ,本 文采 用 章 文 波 "" 提出 的 利用 逐 
月 雨量 估算 降雨 侵蚀 力 的 模型 。 采 用 滑 北 矿区 内 及 
周围 区 域 气象 站 点 的 降雨 数据 ,根据 公式 (2) 计 算 
各 站 点 降雨 侵蚀 力 值 ,并 利用 克 里 金 方法 插值 得 到 


滑 北 矿区 年 降雨 侵蚀 力 栅 格 图 。 


(2) 


式 中 :已 ,为 第 ;年 /月 的 降雨 量 (mm) ;NN 为 年 数 ;R 
为 多 年 平均 降雨 侵蚀 力 (MJ © mm .hm 一 + h^! - 
a`); a B 为 模型 参数 : a =0.183 3,B=1.995 7。 

土壤 可 蚀 性 因子 是 反映 土壤 性 能 和 土壤 被 刨 的 
难 易 程度 的 指标 。 本 文 使 用 在 黄土 高 原 地 区 具有 较 
好 适用 性 的 WILLAMS 等 ”在 侵蚀 力 /生产 力 影响 
模型 EPIC ( Erosion-Productivity Impact Calculator) 中 
发 展 的 天 值 估算 方法 进行 计算 ,公式 如 下 : 


K= fo. 240. 3cl ossi] " 


) x(ro-2 


0.78, 
| 


( 75.51 +22. 98,) 


0. 25C 


5; 
io +S, 


[1.0 - (3) 


S, +e 


式 中 :$。 为 砂粒 (0.1 ~2 mm) 含 量 (% );5; 为 粉 粒 
(0.002 ~0.1 mm) 含 量 (% );C; HA ( <0. 002 


S 
mm) & fi (70) ;C 为 有 机 碳 合 量 (%) 58, =1- 1。 


坡 长 和 坡度 因子 通常 被 称 为 地 形 因 子 (15), 反 
映 了 区 域 地 形 地 貌 特征 对 土壤 侵蚀 造成 的 影响 。 本 
文采 用 在 黄土 高 原 地 区 具有 和 较 好 适用 性 的 Foster 
等 3- 中 建立 的 坡 长 因子 和 坡度 因子 的 计算 方法 , 计 
算 公式 为 : 
L= (x3) 


B= (sing/ 0.089 6) / [3 (sn0) +0.56| (5 


"8/ (1+4) 


10. 8sing +0. 03, 0<5.142 8? 


(4) 


(6) 


S-/(16.8sing -0.5, 5.142 8°<0<14.036 2? (7) 


21.91sin0 -0. 96, 0=14. 036 2? 
式 中 : A 为 坡 长 ;m 为 坡 长 指数 ; 0 为 在 DEM 中 提 
取 的 坡度 记 为 坡 长 因子 ;$ 为 坡度 因子 。 

地 表 植 被 覆盖 因子 (C) 反 映 了 植被 等 地 表 履 盖 
对 土壤 的 保护 作用 ,植被 覆盖 度 因 子 是 指 在 其 他 条 
件 相 同 的 情况 下 ,植被 覆盖 的 土壤 流失 量 与 裸 土 的 
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土壤 流失 量 之 间 的 比值 , 值 在 0 -1 zm] s zx 
采用 蔡 尝 法 等 "建立 的 植被 覆盖 度 与 植被 覆盖 度 
因子 (C) 之 间 的 关系 来 估算 C 值 ,计算 公式 为 : 
l, 0 xf, x0. 196 
C= 0.650 8 -0.343 6lg(f.), 0.196 <f, x78. 396 
0, f, 278.396 


|. NDVI - NDVI,,, 
e ^ NDVI — NDVI; 


(8) 


L 
SSW = Y N,o,,88T = No? (10) 
h=1 


RP :h=1, 4, L AEE Y RAF X 的 分 层 , 即 分 类 
或 分 区 ;N 入 分别 为 屋 h 和 全 区 的 单元 数 ;o? 和 
a^ 分 别 是 层 h 和 全 区 的 了 值 的 方差 。SSW 和 SST 
分 别 为 层 内 方差 之 和 及 全 区 总 方差 。g 为 某 一 自 变 
量 对 区 域 土壤 侵蚀 空间 差异 的 解释 力 , 值 域 为 [0， 
1] ,9 值 越 大 说 明 该 因子 的 作用 越 明显 ,由 此 可 对 研 
究 区 土壤 侵蚀 的 主导 因子 进行 判别 。 

风险 探测 器 用 于 判断 两 个 子 区 域 间 的 属性 均值 


式 中 :人 为 植被 覆盖 度 (% ) ; C 为 地 表 植 被 覆盖 度 
KF ; NDVI , NDVI, f] H] ENVIS. 1 统计 获得 累计 
像 元 值 $% 处 及 95% 处 ,分 别 记 为 裸 土 值 和 植被 完 
全 覆盖 值 。 在 ArcGIS 中 利用 顶 格 计算 器 编写 上 述 
公式 得 到 2018 年 研究 区 地 表 植 被 覆盖 度 因 子 (C ) 
的 栅 格 图 层 。 

水 土 保持 措施 因子 (P) 是 指 采用 水 土 保持 措施 
后 的 土壤 侵蚀 量 与 未 采取 水 土 保持 措施 的 顺 坡 种 植 
时 的 土壤 侵蚀 量 之 比值 ” 。 不 同 水 保 措 施 效果 的 
差异 ,可 以 通过 对 不 同 土地 履 被 类 型 进行 赋值 表达 ， 
但 国内 目前 尚 没有 统一 的 赋值 标准 。 根 据 遥 感 影像 
和 渭 北 矿区 的 实际 情况 ,结合 前 人 的 研究 成 
果 ,对 不 同 土地 利用 类 型 赋予 不 同 的 尸 值 , 确 
定 耕 地 、 林 地 水系、 建设 用 地 、 裸 地 的 P 值 分 别 为 
0.45 .1.0.0.0.2, 得 到 PP 因子 分 布 图 层 。 

2.2.2 地 理 探测 器 ”地理 探测 器 的 原理 是 通过 分 
析 辨 别 要素 不 同 空 间 类 别 分 区 之 间 的 异 质 性 探究 其 
背后 驱动 力 , 在 定量 探测 地 理 现象 的 驱动 因素 方面 
已 经 有 了 一 定 的 应 用 ,并 取得 了 较 好 的 研究 结 
PU- 。 受 到 地 带 和 非 地 带 性 因素 的 影响 ,地 理 要 
素 的 空间 分 布 往 往 并 不 均匀 ,这 种 空间 差异 通过 分 
类 或 者 分 区 进行 表达 称 为 空间 分 层 异 质 性 。 在 地 理 
探测 器 中 通过 层 内 方差 之 和 小 于 层 间 总 方差 对 这 一 
现象 进行 判定 ,分 异性 的 大 小 由 9 值 来 衡量 1。 地 
理 探测 需 的 主要 组 成 部 分 是 因子 探测 器 风险 探测 
器 生态 探测 器 和 交互 作用 探测 器 。 因 子 探测 器 主 
要 用 于 检验 某 种 要 素 是 否 对 因 变 量 的 空间 分 布 格局 
产生 了 影响 ,并 利用 q 值 来 衡量 其 影响 力 的 大 
CU ,其 计算 模型 为 : 


是 否 具 有 显著 差别 ,通过 计算 比较 各 个 影响 因子 在 
不 同类 别 分 区 下 的 土壤 侵蚀 量 均值 来 进行 高 风险 区 
域 的 搜索 ;交互 作用 探测 器 可 以 识别 两 个 自 变 量 因 
子 之 间 的 交互 作用 并 进行 定量 计算 , 即 评估 两 个 影 
响 因 子 共同 作用 时 是 否 会 增加 或 减弱 对 因 变 量 的 解 
释 力 ,通过 计算 比较 单 因 子 q 值 和 交互 两 个 因子 的 
q 值 可 以 分 析 各 因子 的 交互 作用 对 土壤 侵蚀 量 的 影 
响 力 。 相 比 于 其 他 统计 方法 ,这 无 疑 是 地 理 探测 器 
的 最 大 优势 。 

为 了 分 析 渭 北 矿 区 土壤 侵蚀 的 主导 因子 ,利用 
地 理 探测 器 对 植被 覆盖 度 、 多 年 平均 降水 量 、 坡 度 、 
土地 利用 类 型 和 煤炭 年 产量 5 个 因子 进行 分 析 。 由 
于 地 理 探测 器 要 求 输入 的 变量 为 类 别 数据 ,因此 需 
要 对 连续 型 数据 进行 离散 化 。 其 中 土地 利用 数据 根 
据 类 别 进行 分 类 ,结合 王 劲 峰 等 提出 的 数据 离散 
化 方法 及 先 验 知识 "将 多 年 平均 降水 量 数据 按照 
等 间距 分 类 法 分 成 9 类 ,植被 覆盖 度数 据 按照 
<0.3.0.3~0.4.0.4~0.5.0.5~0.6.0.6~0.7、 
0.7 ~0.8.0.8 ~0.9.0.9~1 分 为 8 类 ,坡度 按照 
<5° 5° ~10° 109 ~15°.15° ~20° 909 ~25° 25° ~ 
30° 30° ~35°、>35° 等 级 分 为 8 类 ,煤炭 年 产量 数 
据 在 进行 地 理 位 置 识别 后 手动 添加 ,各 煤矿 点 的 具 
体 数据 如 图 2 所 示 。 利 用 ArcGIS10. 2 将 渭 北 矿区 
划分 为 2 km x2 km 的 研究 单元 ,并 在 每 个 研究 单元 
中 心 点 设置 样本 点 ,将 土壤 侵蚀 量 及 各 影响 因子 分 类 
值 赋 予 创建 的 离散 点 作为 地 理 探测 器 的 运行 数据 。 


3 结果 分 析 


3.1 土壤 侵蚀 评价 

应 用 RUSLE 模型 对 渭 北 矿区 土壤 侵蚀 量 进行 
估算 , 渭 北 矿区 2018 年 土壤 侵蚀 量 的 范围 为 0 ~ 
2 542.46 t hm” «a^! 。 按 照 国家 土壤 侵蚀 强度 分 
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图 2 渭 北 矿区 部 分 煤矿 实地 调研 资料 


Tig.2 Some field research data of coalmines in Weibei mining area 


严重 。 通 过 谷 加 县 级 行政 区 划 可 以 得 到 王 益 区 三 


强度 主要 为 微 度 和 轻 度 侵蚀 ,其 中 微 度 侵蚀 面积 

比 为 83.92% , 轻 度 侵蚀 面积 占 比 为 7.90% 。 中 度 
至 剧烈 土壤 侵蚀 强度 在 消 北 矿区 内 也 有 零星 分 布 ， 
但 面积 占 比 均 较 小 ,分 别 为 4.34% 1. 8596 、1.36% 
和 0.63% 。 由 2018 年 渭 北 矿区 土壤 侵蚀 空间 分 布 
(图 3) 可 以 得 到 , 渭 北 矿区 土壤 侵蚀 具有 空间 异 质 
性 ,侵蚀 较为 严重 的 区 域 位 于 研究 区 西南 部 .中 部 和 
东南 部 ,其 中 尤 以 西南 部 和 中 部 地 区 土壤 侵蚀 最 为 


原 县 北部 .富平 县 北部 .耀州 区 东部 .白水 县 浦城 县 
北部 OREL VS 合 阳 县 北部 和 韩 城市 中 部 地 区 的 
土壤 侵蚀 较为 严重 。 对 比 区 域 地 形 图 可 以 发 现 ,这 
些 地 区 地 势 起 伏 较 大 ,地 形 较 为 复杂 , 且 有 清河 . 石 
川 河 白水 河北 洛 河 等 河流 经 过 。 由 于 煤炭 开采 等 
原因 降低 了 这 些 地 区 的 植被 覆盖 度 ,使 得 植物 对 土 
坏 的 保护 作用 减弱 ,呈现 出 较为 严重 的 土壤 侵蚀 现 
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Fig.3 Spatial distribution of soil erosion in the Weibei mining area in 2018 
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3.2 基于 地 理 探测 器 的 土壤 侵蚀 定量 归 因 

3.2.1 土壤 侵蚀 影响 因子 显著 性 研究 ”因子 探测 
器 运行 结果 表明 ,不 同 影响 因子 对 土壤 侵蚀 量 的 解 
释 力 具有 显著 差异 。 从 渭 北 矿区 全 范围 来 看 ,各 影 
响 因 子 对 土壤 侵蚀 量 的 解释 力 从 大 到 小 依次 为 植被 
覆盖 度 多 年 平均 降水 量 、 土 地 利用 类 型 坡度 和 煤 
炭 年 产量 。 表 1 为 对 所 有 因子 的 9 值 计算 结 果 , 结 
果 表 明 植 被 覆盖 度 因 子 具 有 最 高 的 4 值 ,说 明 在 这 
些 影响 因子 中 植被 覆盖 度 是 决定 渭 北 矿区 土壤 侵 伺 
空间 格局 的 主导 因子 。 煤 炭 年 产量 因子 的 解释 力 最 
小 , 且 未 通过 显著 性 检验 ,说 明 该 因子 对 土壤 侵蚀 的 
影响 程度 最 低 。 根 据 实地 调研 的 结果 , 渭 北 矿区 作 
为 陕西 省 最 早 建 成 的 煤炭 生产 基地 ,经 过 长 期 开采 ， 
矿区 已 经 出 现 资源 枯竭 等 诸多 问题 ,多 数 煤矿 已 经 
闲置 或 关闭 ,煤矿 数量 大 幅 减 少 日 实施 了 复 绿 工程 。 
此 外 渭 北 矿区 煤矿 多 为 井 工 开采 方式 , 相 比 于 露天 
煤矿 对 土壤 的 扰动 程度 降低 ,因此 通过 分 析 2018 年 
的 数据 发 现 ,煤炭 开采 已 不 是 渭 北 矿区 土壤 侵蚀 的 
主要 影响 因子 。 对 比 自然 环境 因素 与 人 为 因素 的 解 
释 力 可 以 发 现 , 植 被 覆盖 度 .多 年 平均 降水 量 等 自然 
因素 的 q 值 大 于 煤炭 年 产量 等 人 为 因素 的 g 值 ,说 
明 自 然 环境 因素 是 影响 渭 北 矿区 土壤 侵蚀 空间 分 布 
的 主要 驱动 力 。 

3.2.2 土壤 侵蚀 高 风险 区 域 识 别 ”风险 探测 器 结 
果 如 表 2 所 示 , 土 地 利用 类 型 中 裸 地 发 生 土壤 侵蚀 
的 风险 最 高 ,这 是 因为 渭 北 矿区 地 处 黄土 分 布 区 域 ， 
在 水 力 等 作用 下 极 易 发 生 侵 蚀 。 裸 地 环境 条 件 较 
差 , 少 有 植物 生长 ,同时 科学 有 效 的 水 土管 理 措施 实 
施 并 不 到 位 ,造成 士 壤 侵蚀 现象 较为 严峻 。 通 过 对 
多 年 平均 降水 量 分 层 土壤 侵蚀 均值 进行 比较 发 现 降 
水 量 在 536. 42 ~ 548. 62 mm 间 的 区 域 发 生 土 壤 侵 
刨 的 风险 最 高 。 又 加 渭 北 矿区 多 年 降水 量 空间 分 布 
图 可 以 得 到 ,降水 量 在 该 区 间 范 围 的 区 域 主要 位 于 
研究 区 中 部 (白水 县 Oe EL i e E) ,这 些 地 区 的 
土壤 侵蚀 现象 均 较为 严重 。 通 过 对 比 不 同 坡度 分 区 
对 应 的 土壤 侵蚀 量 均 值 可 以 发 现 随 着 坡度 的 增加 ， 
土壤 侵蚀 量 呈 先 增加 后 减少 的 趋势 ,在 坡度 20? ~ 
25" 区 间 出 现 抛 点 , 即 为 表 2 中 坡度 的 高 风险 区 。 在 

表 1 各 影响 因子 4 值 统计 


Tab.1 q value of influencing factors 


解释 力 土地 利 ”多 年 平均 降 ”坡度 植被 覆 煤炭 年 产量 
T^^ 用 类 型 ”水 量 /mm —/ 3H 7/10 1:2 
qti 0.00172 — 0.0594 0.0169 0.1142 0.0103 


R2 各 影响 因子 土壤 侵蚀 高 风险 区 域 及 其 
平均 值 /t. hm? . a^! 


Tab.2 High risk areas of soil erosion and its mean 


-2 -1 


value in influencing factors / t- hn » a 
土地 多 年 平 ea ”植被 煤炭 年 
县 坡度 E- ER 

利 降水 量 y T us 产量 

类 型 / mm RE /10*t-a^! 
高 风险 区 b 536.42 -548.62 20-25 «0.3 62 
土壤 侵蚀 13.03 19. 22 13.77 29.32 一 
均值 
注 :“ - ”表示 未 通过 置信 水 平 95% 的 检验 


对 植被 覆盖 度 进行 层 间 分 析 后 发 现 , 渭 北 矿区 内 植 
被 覆盖 度 与 土壤 侵蚀 量 密切 相关 ,土壤 侵蚀 量 随 着 
植被 覆盖 度 的 升 高 而 减少 ,植被 覆盖 度 小 于 0.3 时 
土壤 侵蚀 量 均 值 达到 最 大 值 。 又 加 植被 覆盖 与 土壤 
侵蚀 空间 分 布 图 发 现 , 侵蚀 程度 较 轻 的 区 域 主要 分 
布 在 高 植被 覆盖 区 ,侵蚀 程度 较 严 重 的 区 域 主要 分 
布 在 低 植被 覆盖 区 ,可 见 植 被 对 土壤 侵蚀 的 抑制 作 
用 较为 显著 。 

3.2.3 土壤 侵蚀 影响 因子 交互 作用 研究 ”交互 作 
用 探测 器 统计 结果 如 表 3 所 示 , 可 见 两 种 影响 因子 
的 协同 作用 会 增强 对 土壤 侵蚀 的 解释 力 。 对 统计 结 
果 进 行 分 析 可 以 得 到 , 除 降 水 与 坡度 的 交互 作用 外 ， 
交互 作用 解释 力 排 在 第 一 位 的 均 为 植被 覆盖 度 与 其 
他 因子 的 协同 作用 ,是 土壤 侵蚀 的 显著 控制 因子 ,这 
意味 着 植被 覆盖 度 相差 较 大 的 区 域 往往 土壤 侵 刨 程 
度 有 着 明显 差异 。 将 各 因子 与 植被 覆盖 度 交 互 作用 
后 的 g 值 与 单 因子 q 值 进行 比较 ,坡度 与 植被 覆盖 
度 因 子 交 互 作用 后 q 值 变化 幅度 最 大 ,二 者 协同 作 
用 后 对 土壤 侵蚀 的 解释 力 约 是 坡度 单 因 子 解 释 力 的 
25 倍 。 因 此 应 当 注 重 坡 面 恢复 植被 ,植树 种 草 , 增 
加 入 渗 , 减 少 地 表 径 流 和 土壤 的 破坏 转移 。 


4 结论 


通过 对 渭 北 矿区 2018 年 土壤 侵蚀 量 进行 估算 


表 3 土壤 侵蚀 影响 因子 交互 作用 9 值 统 计 
Tab.3 q values of dominant interactions between 


two covariates 


2 Bb AE EI 
土地 利 0.017 0.082 0.076 0.240 0.039 
降水 0.082 0.059 0.226 0.215 0.082 
坡度 0.076 0.226 0.017 0.435 0.036 


植被 覆盖 度 0.240 0.215 0.435 0.114 0.138 
_ 煤炭 年 产量 — 0.099 — 0.080 — 0.036 0.138 — 0.020 _ 
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以 及 对 驱动 因素 进行 分 析 后 得 出 以 下 结论 : 

(1) 滑 北 矿区 土壤 侵蚀 强度 主要 为 微 度 和 轻 度 
侵蚀 ,侵蚀 严重 的 地 区 主要 位 于 研究 区 西南 部 .中 部 
和 东南 部 。 在 县 级 行政 区 单元 来 看 ,侵蚀 较为 严重 
的 地 区 主要 在 耀州 区 东部 富平 县 西部 .白水 县 东 
部 泪 城 县 西部 、 合 阳 县 南部 和 韩 城市 中 部 。 通 过 县 
加 河流 和 植被 覆盖 分 布 图 可 以 发 现 ,这 些 区 域 大 多 
有 河流 经 过 且 植 被 覆盖 度 相 对 较 低 。 

(2) 通过 利用 地 理 探测 器 对 土壤 侵蚀 各 影响 因 
子 进行 分 析 发 现 植被 覆盖 度 的 解释 力 最 大 ,说 明 植 
被 覆盖 度 是 影响 渭 北 矿区 土壤 侵蚀 空间 分 布 的 主导 
因子 。 在 对 各 因子 进行 层 间 侵蚀 量 差 异 进行 分 析 后 
识别 出 渭 北 矿区 土壤 侵蚀 高 风险 区 位 于 坡度 20° ~ 
25° .植被 覆盖 度 < 0.3 和 用 地 类 型 为 裸 地 的 区 域 ， 
因此 这 些 区 域 需 注重 实施 水 土 保持 措施 。 

(3) 由 于 土壤 侵蚀 是 多 个 因子 共同 作用 的 结 
TR ,进一步 对 各 因子 进行 交互 作用 研究 发 现 各 影响 


因子 两 两 交互 作用 对 土壤 侵蚀 的 解释 力 均 大 于 单 因 
子 解释 力 , 且 除去 植被 覆盖 度 与 降水 的 交互 作用 解 
释 力 略 小 于 坡度 与 降水 的 交互 作用 解释 力 外 ,植被 
窗 盖 度 因 子 与 其 他 因子 的 协同 作用 解释 力 均 排 在 第 
一 位 ,是 渭 北 矿区 土壤 侵蚀 的 显著 影响 因子 。 

(4) 渭 北 矿区 地 处 黄土 高 原 ,特殊 的 自然 环境 
和 强烈 的 人 类 活动 都 加 重 了 区 域 土壤 侵蚀 。 在 今后 
的 水 土 保持 工作 实施 过 程 中 除 进一步 对 煤矿 区 进行 
生态 修复 外 ,还 应 注重 坡 面 恢复 植被 ,进而 改良 土壤 
增加 入 渗 , 同 时 注重 开展 对 地 形 复 杂 地 区 的 水 土 保 
持 工 作 。 受 到 收集 到 的 数据 的 精度 限制 ,所 选用 的 
遥感 影像 分 辨 率 为 30 m x30 m, 加 之 研究 区 域 面积 
较 小 ,因此 在 后 续 的 研究 中 应 当 着 重 提高 数据 源 的 
精度 ,进而 得 到 更 好 的 研究 成 果 。 
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Evaluation of soil erosion and driving factors analysis in Weibei mining area 


ZOU Ya-jng, YAN Qing-wu, TAN Xue-ling, WANG Jin 
( College of Environment Science and Spatial Informatics , China. University of Mining and Technology , Xuzhou ,221116 , Jiangsu , China) 


Abstract: The Weibei mining area in Shaanxi Province , China is located in the Loess Plateau. Affected by natural 
environment and human activities such as coal mining, soil erosion in this area is more serious and the ecological 
environment is more fragile. In this paper,the simulation of soil erosion and its quantitative attribution analysis have 
been conducted in the mining area based on the RUSLE model and the geographical detector method. The influen- 
cing factors ,such as vegetation cover, perennial average precipitation, slope, land use type and annual coal output 
have been taken into consideration. The results of this study provided the advice on the soil and water conservation 
in mining areas. The results indicate as follows; (1) The intensity of soil erosion in general was mainly slight and 
mild in Weibei mining area. The soil erosion was severe in the central , southeastern and southwestern regions of the 
study area. (2) Vegetation coverage and perennial average precipitation are the dominant factors of soil erosion in 
Weibei mining area. Moreover ,the conditions that are likely to induce soil erosion are; slopes between 20? and 25? , 
vegetation coverage less than 0.3 and bare lands. (3) The synergistic effect of the factors on soil erosion is greater 
than that of single factor. Therefore ,the combined action of multiple factors have a significant impact on soil erosion 
in Weibei mining area. 
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